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As disc galaxies form in a hierarchical cosmology, massive merging satellites are preferentially dragged towards the disc plane. 
The material accreted from these satellites forms a dark matter disc that contributes 0.25 - 1.5 times the non-rotating halo density 
at the solar position. Here, we show the importance of the dark disc for indirect dark matter detection in neutrino telescopes. 
Previous predictions of the neutrino flux from WIMP annihilation in the Earth and the Sun have assumed that Galactic dark matter 
is spherically distributed with a Gaussian velocity distribution, the standard halo model. Although the dark disc has a local density 
comparable to the dark halo, its higher phase space density at low velocities greatly enhances capture rates in the Sun and Earth. 
For typical dark disc properties, the resulting muon flux from the Earth is increased by three orders of magnitude over the SHM, 
while for the Sun the increase is an order of magnitude. This significantly increases the sensitivity of neutrino telescopes to fix 
or constrain parameters in WIMP models. The flux from the Earth is extremely sensitive to the detailed properties of the dark 
disc, while the flux from the Sun is more robust. The enhancement of the muon flux from the dark disc puts the search for WIMP 
annihilation in the Earth on the same level as the Sun for WIMP masses < 100 GeV. 
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1. Introduction 

A mysterious dark matter makes up most of the mass in 
the Universe, providing us with a unique window into physics 
beyond the Standard Model. Among the many plausible dark 
matter particle candidates, Weakly Interacting Massive Parti- 
cles (WIMPs) CI S S that arise in theories with supersym- 
metry (SUSY) [4] or universal extra dimensions (UED) |5j,|6[], 
stand out as well-motivated and detectable. WIMPs may be 
detected directly by scattering in laboratory detectors J3l, or in- 
directly by their annihilation products. The annihilation rate 
scales as the square of the WIMP density, so the most lumi- 
nous sources are expected to be near the Galactic centre or the 
centre of dark matter sub-halos, where the dark matter density 
peaks |1 11 11 [13, [ll El]. In addition, the Sun and Earth cap- 
ture WIMPs and may be seen as sources of WIMP annihilation 
lfl3i FA [l5l Il6l rTvh . In all cases, the annihilation rate is sensi- 
tive to the dark matter's phase space structure. 

WIMPs can be gravitationally trapped inside the Sun and 
Earth by elastic scattering, if the final WIMP states have ve- 
locities below the escape velocity. To date, annihilation rates 
in the Sun and Earth have been estimated using the Standard 
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Halo Model (SHM), which is modelled as a smooth, spherically 
symmetric density component with a non-rotating Gaussian ve- 
locity distribution (sometimes with an anisotropic velocity dis- 
persion tensor). However, recent dark matter-only simulations 
show deviations from this picture in the form of many small 
amplitude density fluctuations on ~ 100 pc scales, a few large 
amplitude density fluctuations due to the presence of substruc- 
ture (which makes up ~ 0.5% of the mass at the solar circle), 
and relatively small deviations of the velocity distribution from 
Gaussian 1191 12011 . Even if such structures were to survive in 
the presence of a baryonic disc and bulge, the indirect detec- 
tion signal from the Earth and Sun is unlikely to deviate signif- 
icantly from the SHM prediction. This is because the annihi- 
lation rate is sensitive to the phase space density averaged over 
long (> 100 Myr) timescales. As a result, indirect detection by 
annihilation in the Earth and Sun is only sensitive to the local 
dark matter macro structure. 

There is at least one local macro-structural component be- 
yond the SHM, which has been discovered in recent simulations 
of galaxy formation that include baryons. The baryonic disc of 
the Milky Way draws satellites closer to the disc plane by dy- 
namical friction, where they are disrupted by tides [21]. This 
results in a thick disc of dark matter with a mid-plane density 
of 0.25-1.5 times the local dark halo density and kinematics 
similar to the thick disc of stars 11221 p23n . The dark disc boosts 
the flux in direct detection experiments at low energies and in- 
creases the annual modulation signal with an energy-dependent 
pha se shift that will betray the mass of the dark matter particle 
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In this Letter, we focus on the importance of the dark disc 
for the detection of neutrinos from WIMP annihilation in the 
cores of the Sun and Earth. WIMPs can annihilate into a wide 
range of final products, of which muon neutrinos can escape 
and reach terrestrial detectors. On Earth, these muon neutri- 
nos produce muons in charged current interactions with nuclei: 
+ N — » fT + X. The ultra-relativistic muons can be detected 
by their Cerenkov radiation in large water or ice-based neu- 
trino telescopes. So far, neutrino telescopes have found no ev- 
idence for high-energy neutrinos of astrophysical origin above 
the detected atmospheric neutrino background. The most strin- 
gent bounds on high-energy neutrinos from the Sun and Earth 
come from Super-Kamiokande l25ll . AMANDA i26ll and Ice- 
Cube [27]. Super-Kamiokande, a 50 kt water Cerenkov detec- 
tor, and AMANDA, located in the ice sheet at the Amundsen- 
Scott South pole station, have been taking data throughout the 
past decade. Their muon flux limits from the Earth and the Sun 
are of order O^, ~ 10 3 -10 4 kirT 2 yr~' for energies E M > 1 GeV, 
where the tighter constraints apply to higher WIMP masses. 
These flux limits and those shown in the figures below are based 
on the assumption of a hard WIMP annihilation spectrum. Ice- 
Cube, currently under construction at the site of the AMANDA 
experiment, has achieved a similar sensitivity with even a small 
fraction (~ 27%) of the construction completed, and is expected 
to have a 5-year sensitivity to flux from the Sun of ~ 70 
km 2 yr~' above a WIMP mass of ~ 200 GeV. The expected 
5 -year sensitivity to flux from the Earth is (D^ ~ 20 krrT 2 yr~' . 

In this Letter, we show that the presence of a dark disc dra- 
matically increases the parameter space of WIMP models to 
which neutrino telescopes are sensitive. In $2] we introduce the 
dark disc and its properties used in this Letter. In §|3] we de- 
rive the capture and annihilation rates of WIMPs in the Sun 
and Earth. In §@] we describe a particular WIMP candidate 
- the neutralino in the Constrained Minimal Supersymmetric 
Model (CMSSM). In §0 we calculate the expected muon flux 
for CMSSM neutralinos that are consistent with both astrophys- 
ical and collider constraints, and we show how the flux changes 
once we include the dark disc. In ^6] we discuss how our re- 
sults are affected by uncertainties in the dark disc phase space 
density. Finally, in $7] we present our conclusions. 

2. The dark disc 

As disc galaxies form in a hierarchical cosmology, many 
smaller satellite galaxies are accreted and dissolve. Read et al. 
1 22] looked at the interaction of a disc with a cosmological dis- 
tribution of satellites. They found that low inclination massive 
satellites were preferentially dragged into the disc plane by dy- 
namical friction. The material accreted from these satellites set- 
tled into a thick disc of dark matter. Read et al. lEill studied ab 
initio simulations of galaxy formation where a thick dark disc 
was also found. In both studies, the dark disc had a density 
Pd that is ~ 0.25 - 1.5 times the local density of dark matter 
in the SHM ph, a rotation lag with respect to the stellar disc 
of |v Q | ~ — 150km/s, and a near-isotropic ID dispersion of 
cr ~ 50 - 90km/s. The uncertainties in these numbers reflect 
the unknown stochastic merger history of our Galaxy. 



An accreted thick stellar disc forms concurrent with the 
dark disc and shares similar kinematics. Although the origin 
of the Milky Way thick disc of stars remains under investiga- 
tion, it has kinematic properties remarkably similar to the me- 
dian dark disc distribution with | v Q | =40-50 km/s with respect 
to the local circular velocity, and dispersions of (<x s , cr^, cr z ) = 
(63,39,39)km/s JUS. For our present study, we model the 
distribution functions of both the SHM and dark disc as Gaus- 
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where u is the heliocentric WIMP speed, cr is the ID velocity 
dispersion and v© is the lag between the dark matter particles 
and the Sun. p is the WIMP density at the solar circle and M x 
is the WIMP mass. For the SHM, |v | = 220 km/s and cr = 
|v |/ V2. For the dark disc, we assume fiducial values of: pa = 
Ph and |v Q | = cr — 50 km/s, consistent both with the median of 
the ranges found in [22] and [23], and with the kinematics of 
the Milky Way thick disc stars. The sensitivity of our results to 
these parameters is discussed in §|6] 

3. Capture from the dark disc 

The capture rate from a nuclear species / per unit volume 
shell of a celestial body is given by [ 14] 
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where /(u) is the velocity distribution normalised such that 
J f(u)d 3 u = p/M x . The velocity w at a given shell is related to 
the velocity at infinity u and the escape velocity v at the shell 
by w — Vm 2 + v 2 . The WIMP-nucleus cross section is <x,-, and 
yii is the number density of nuclear species i. The upper limit of 
the integration is 
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where m-, is the mass of nuclear species 2. This ensures that 
only WIMPs that can scatter to a velocity lower than the escape 
velocity v are included. For fixed mass m, and escape velocity 
v this upper cut off decreases with increasing M x . 

The annihilation rate per unit volume of WIMPs in the body 
is given by 

^ =<cr A v) n 2 (f,x), (4) 

where (oa v)o is the velocity-averaged annihilation cross section 
in the limit of non-relativistic speeds, and n x (t, x) is the number 
density of WIMPs in the body. If WIMPs quickly thermalize 
with nuclei in the body once captured, the number density of 
WIMPs in that body can be described by 

n x {t,x) = N x (t)n x (x), (5) 

where f dVn x (x) = 1. In that case, the number N x (t) of WIMPs 
in the body is given by the solution to 



N x (t) = C-2T A 



(6) 
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if M x » ;«,, where the total capture rate is C — 2, C,- IU8U15I1 . 



The factor of 2 in Eq. [6]reflects the fact that for self-annihilating 
particles, two WIMPs are destroyed in each annihilation. If the 
total capture rate C is constant with time, the annihilation rate 
is given by 



r A = | tanh 2 (f/r), 
with the equilibrium time r given by 



r = (CC A y 1/2 , C A = 2<(rv>o I dVnt(x) 
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For a WIMP of M^, ~ 100 GeV with purely spin-independent 
interactions with baryons, C is about nine orders of magnitude 
greater for the Sun than for the Earth, while Ca is about three 
orders of magnitude smaller [15]. Thus, equilibrium timescales 
tend to be orders of magnitude shorter in the Sun than in the 
Earth. If the spin-independent WIMP-proton cross section is 
of = 10~ 43 cm 2 , t ~ 10 8 yr for the Sun and t ~ 10 11 yr 
for the Earth [4]. Therefore, if the age of the Solar system 
is t Q ~ 4.5 Gyr, ? q /t » 1 in the Sun and F A - C/2 and is 
constant. For this set of WIMP parameters, WIMP annihila- 
tion will have reached equilibrium in the Sun. However, in the 
Earth, t /r «: 1, and F^ oc C 2 and is growing with time. For 
some of the models that we consider, the interaction rate be- 
comes low enough that equilibrium is broken in the Sun also. 
However, these models appear near the bottoms of the plots, far 
away from the interesting detection thresholds. 

4. The dark matter candidate 

For detailed calculations of capture and annihilation rates, 
we must assume a specific WIMP model for the particle's mass, 
scattering cross sections and annihilation channels. In super- 
symmetric extensions of the Standard Model with conserved R- 
parity, the lightest supersymmetric particle is a natural WIMP 
candidate. We choose the lightest neutralino in the CMSSM 
as the dark matter particle. The CMSSM reduces the free pa- 
rameters in supersymmetry to three parameters at a gauge uni- 
fication scale: the gaugino mass ;«i/2, the scalar mass mo and 
the tri-linear coupling Aq. At the electroweak scale, the ratio 
of the Higgs vacuum expectation values tan(y3) and the sign of 
the Higgs/higgsino mass parameter yi are selected. Only the 
square of p is calculated from the potential minimisation con- 
ditions of electroweak symmetry breaking. We scanned the 
CMSSM parameter space in the range of mi/2 e [0 4]TeV, 
m € [0 4] TeV, A € [-7 7] TeV and tan(/3) e [20 65] for 
H > 0. The allowed parameter space, consistent with current 
experimental constraints is found with the publicly available 
SuperBayes package 12811 . This Markov Chain Monte Carlo 
(MCMC) algorithm calculates the Bayesian posterior probabil- 
ity at each parameter point from the compatibility of the theo- 
retical predictions with experimental constraints. The most rel- 
evant experimental constraints used in the MCMC are the cos- 
mologically allowed relic density measured by the Wilkinson 
Microwave Anisotropy Probe (WMAP), electroweak precision 
observables and limits on the Higgs and lightest neutralino mass 



from colliders. All constraints listed in Ruiz de Austri, Trotta & 
Roszkowski 12811 have been used, except that the value for the 
dark matter density has been updated to the WMAP 5 year data 
release value, Cl DM h 2 = 0.1099 ± 0.0062, where h is the Hub- 
ble constant in units of 100/;km/s/Mpc and h = 0.719 is their 
best fit value 1E9I1 . Flat priors are used in the calculation of the 
Bayesian probability. Although the resulting allowed parame- 
ter space depends on our choice of priors, we are interested in 
the relative change in flux caused by the addition of the dark 
matter disc and thus our choice of priors is not important. The 
chains used in this study contain a total of 0.9 x 10 6 samples. 
The posterior probability density functions presented below are 
normalised to their maximum values, and not to a total proba- 
bility of 1. Accordingly we label these as "relative probability 
densities". An interface of the SuperBayes MCMC algorithm to 
DarkSusy (v5.03) Q routines is used to calculate the capture 
rate and the resulting muon flux at the Earth. 

5. Neutrino flux from the Sun and Earth 

5.7. The Sun 

In Fig. 1 we show the muon flux with an energy above 
1 GeV at the Earth resulting from WIMP capture and annihi- 
lation in the Sun. We show the flux as a function of M x for 
the SHM (left) and the dark disc (right) for pa/ph = 1 and 
Ph - 0.3 GeV/cm 3 . The higher phase space density at low 
velocities for the dark disc strongly enhances the capture rate 
and hence the resulting muon flux at the detector site. The flux 
expected from the dark disc is larger by approximately an or- 
der of magnitude (depending on the specific model) compared 
to the flux expected from the SHM, since the capture rate in- 
creases by approximately an order of magnitude, and f /r » 1 
for most of the models in the figure (see $3). 

5.2. The Earth 

Finding the muon flux from WIMP annihilation in the Earth 
is somewhat more complicated. Because the escape velocity of 
the Earth is small (v m 15km/s at the centre), capture is only 
possible for low speed WIMPs unless the WIMP mass is nearly 
identical to that of one of the nuclear species in the Earth [Eq. 
[3j see also LL4j,l32[]. Moreover, the capture rate is disproportion- 
ately sensitive to the lowest speed WIMPs since those WIMPs 
may be captured anywhere in the body, whereas higher speed 
WIMPs may only be captured at the centre where the escape 
velocity is largest. However, the low speed tail of the WIMP 
speed distribution is not precisely known; for speeds relative to 
the Earth of u < 72km/s (the speed of a WIMP at the escape 
velocity from the Solar system, moving in the direction oppo- 
site to the Earth), the phase space may be occupied by WIMPs 
bound to the Solar system as well as Galactic WIMPs stream- 
ing through the Solar system on unbound orbits. Thus, the an- 
nihilation rate of WIMPs in the Earth depends on the density 
of WIMPs bound to the Solar system, which has not yet been 
definitively determined. 
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Figure 1 : Muon flux <D^ for £ ;1 > 1 GeV at the Earth's surface as a function of M x from neutrinos originating in the Sun, for the SHM (left panel) and the dark disc 
(right panel). The dark disc boosts the muon flux by an order of magnitude for pj/pi, = 1. Current experimental constraints on the muon flux from the Sun from 
Super-Kamiokande [25 ] , AMANDA-II (25l3lh and IceCube22 l27ll along with the expected sensitivity of IceCube80 are over-plotted on the left panel. The closed 
contours show - 95% (red/dashed) and 68% (green/solid) - of the probability density of CMSSM models consistent with both astrophysical and collider constraints, 
and assuming flat priors. The colour-bar gives the relative probability density (see ^4]for details). 




Figure 2: The phase space density at low velocities for the SHM (left panel) and the dark disc (right panel). The black/solid curve is the summed distribution of 
bound (red/dotted) and unbound (orange/dash-dotted) particles from the Solar system simulations, used in the calculation of the capture rates. The blue/dashed line 
shows the distribution of the free space Gaussian approximation. Note the vertical scales of the two plots differ by two orders of magnitude. 



5.3. Sensitivity to the population bound to the Solar system 

There is a spread in predicted bound WIMP distributions. 
Although the following studies predicted the bound WIMP dis- 
tribution for the SHM, the results generalise to arbitrary dark 
matter distributions. On the high end, Gould 113 311 argued that 
the low speed WIMP distribution resulting from gravitational 
capture of WIMPs by the planets should be approximately the 
free space Gaussian distribution function of Eq. ([TJ. This ar- 
gument was based on treating WIMP-planet encounters as lo- 
cal, with the cumulative changes to WIMP speed treated in the 
random walk approximation. Also using the local approxima- 
tion, Lundberg and Edsjo [ 34] found a smaller low speed WIMP 



distribution if they treated the Sun as being infinitely optically 
thick to WIMPs. Damour and Krauss [35] considered a pop- 
ulation of long-lived WIMPs captured in the Solar system by 
elastic scattering in the Sun, but neglected subsequent scatters 
of those WIMPs with solar nuclei. Bergstrom et al. l36Tl found 
that this population could boost the annihilation rate of WIMPs 
in the Earth by a factor of ~ 100 for 60 GeV < M x < 130 GeV. 



More recently, Peter l37l|38D has simulated ~ 10 10 WIMPs 
bound to the Solar system by either gravitational capture or 
elastic scattering in the Sun. Orbits were integrated in a toy 
Solar system consisting of the Sun and Jupiter. WIMP trajecto- 
ries were followed using a modified symplectic integration al- 
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Figure 3: Muon flux for £„ > 1 GeV at the Earth's surface as a function of M x from neutrinos originating in the Earth. The top row is calculated with the free 
space Gaussian approximation of the velocity distribution, while for the bottom row the velocity distribution from the Solar system simulations is used. The sharp 
decrease in the muon flux at high WIMP masses in the bottom row is caused by the kinematic cut-off of the capture rate (see text). Compared to flux from the 
SHM (left panels) the flux from the dark disc (right panels) is boosted by two to three orders of magnitude for pd/Ph = 1, depending on the specific model. Current 
experimental constraints on the muon flux from the Earth from Super-Kamiokande | 25] and AMANDA-II [26, 31] along with the expected sensitivity of IceCube80 
are compared to the flux expected from the SHM. The closed contours show - 95% (red/dashed) and 68% (green/solid) - of the probability density of CMSSM 
models consistent with both astrophysical and collider constraints, and assuming flat priors. The colour-bar gives the relative probability density (see jPlfor details). 
Note the vertical scales of the two bottom plots differ by two orders of magnitude as compared with the top plots. 



gorithm, allowing for the possibility of further elastic scattering 
in the Sun. The orbits were integrated until the WIMPs were 
ejected or scattered onto orbits that no longer intersected the 
Earth. The phase space density distribution of bound WIMPs 
as a function of speed relative to the Earth is shown in Fig. 2 for 
both the SHM and the dark disc. Also shown in Fig. 2 are the 
phase space density distributions of only the Galactic WIMPs 
(unbound to the Solar system) and the free space Gaussian ap- 
proximation (denoted as "Free space" on relevant figures). 

The bound velocity distribution is significantly smaller than 
predicted by Gould [33] and Damour and Krauss 113 511 . and sim- 
ilar to that found by Lundberg and Edsjo yjj] . While part of 
the difference is due to elastic scattering in the Sun (especially 
for the Damour and Krauss population), part of the difference 



is due to simulating orbits in a toy Solar system. The cut-off in 
the velocity distribution at u ~ 9 km/s owes to the phase space 
below being inaccessible to WIMPs in the toy Solar system due 
to the conservation of the Jacobi integral of motion; interaction 
with the inner planets is required to populate lower speeds. This 
cut-off in speed translates to a cut-off in muon flux above a par- 
ticular WIMP mass. Solving Eq. (f3]l for M x , and setting m, to 
the mass of 56 Fe (the dominant atomic species in the core of the 
Earth), we find that the muon flux is exactly zero for M x > 700 
GeV. The impact of the WIMPs bound to the Solar system on 
the WIMP annihilation rate in the Earth depends crucially on 
how effective the inner planets are at populating the phase space 
below u = 9 km/s. 

Given the uncertainty in the low speed WIMP distribution, 
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Figure 4: Total muon flux for > i GeV at the Earth's surface as a function of M x from neutrinos originating in the Sun (left panel) and Earth (right panel). 
In both cases the muon flux is dominated by the dark disc component. Current experimental constraints on the muon flux from the Earth and Sun from Super- 
Kamiokande [25], AMANDA-II l26l[3lll and IceCube22 l27ll are shown. The enhanced flux allows these experiments to constrain a much larger portion of the 
CMSSM parameter space. The projected sensitivity of the IceCube80 experiment will probe a significant fraction of the allowed parameter space. The closed 
contours show - 95% (red/dashed) and 68% (green/solid) - of the probability density of CMSSM models consistent with both astrophysical and collider constraints, 
assuming flat priors. The colour-bar gives the relative probability density (see ^4]for details). 



we calculate annihilation rates for both the Peter 1371 13811 distri- 
bution functions and the free space Gaussian distribution func- 
tion. These span a likely range of the true distribution func- 
tion of low speed WIMPs. The Damour and Rrauss [35] solar- 
captured distribution function is larger than the Gaussian distri- 
bution function for 30 km/s < u < 50 km/s. However Peter 13711 
finds that subsequent scattering in the Sun reduces the lifetime 
and phase space density of these WIMPs below the Gaussian 
phase space density. 

In Fig. 3, we show the muon fluxes from WIMP annihila- 
tion in the Earth for both the SHM and the dark disc assuming 
a muon energy threshold of 1 GeV. For both distribution func- 
tions, the flux from the dark disc is two to three orders of mag- 
nitude above the SHM if pd/ph = 1. This large increase is due 
to the fact that f /r < 1, such that T oc C 2 . Thus, an increase in 
the capture rate of WIMPs in the Earth has a more dramatic ef- 
fect on the muon flux than a similar enhancement in the capture 
rate of WIMPs in the Sun. 

The predicted flux from WIMPs with M x > 100 GeV is 
quite sensitive to the low speed phase space density distribu- 
tion. For the distribution function from the Solar system simu- 
lations, we find the steep drop in flux due to the kinematic cut- 
off in the capture rate for M x > 500 GeV. As a consequence, 
while the enhancement of the muon flux from the dark disc 
puts the search for WIMP annihilation in the Earth on the same 
level as the Sun for M x < 100 GeV, the prospects for detecting 
WIMPs of higher masses is unclear. Precision estimates of the 
low speed tail of the WIMP velocity distribution are necessary 
to determine the prospects for high mass WIMPs. 



6. Discussion 

In Fig. 4, we show the total flux from the Sun and the 
Earth (including capture from both SHM and dark disc com- 
ponents assuming pa/ph = 1) along with current experimental 
constraints. The flux in both cases is dominated by the dark 
disc component. To be conservative, we show the lower bound 
of the expected muon flux from the Earth obtained using the 
phase space density distribution from the Solar system simula- 
tions. The inclusion of the dark disc component significantly 
improves the constraints on the allowed parameter space from 
current experiments. Large area neutrino telescopes such as 
IceCube will be sensitive to a large fraction of the allowed pa- 
rameter space, providing a complementary search for dark mat- 
ter to direct detection experiments. 

Systematic uncertainties owing to the unknown density and 
velocity distribution of the dark disc are especially large for 
the Earth owing to high powers of these parameters in the cal- 
culation of the annihilation flux. For the results presented in 
this Letter, we used pd/ph = 1 and <x rf = 50km/s with the 
mean lag |v Q | = <Jd- For the Earth the dependency is given by 
exp(-|v | 2 /cr 2 )(pd/ph) 2 /cr^ for masses M x > lOOGeV, since 
the part of the WIMP phase space relevant for capture scales as 
exp(-|v | 2 /(2cr 2 ))(p d /ph)/cr^ and the flux depends on the cap- 
ture rate squared. For the Sun the actual dependency is more 
complex 0. An empirical estimate for |v Q | = cr^ of the depen- 



The capture rate scales as 
(Pd/Ph) 
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(9) 



where v cut ~ 2000 ^JGeV/M x km/s is approximately the maximum speed of 
WIMPs which can be captured in the Sun. 
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dency gives a scaling of (pd/Ph)/o"^ with a e [12] depending on 
the particle's mass. This scaling differs from the scaling in the 
Earth owing to the flux being proportional to one power of the 
capture rate. Simulations have shown that all disc galaxies will 
have a dark disc, but the cosmic variance in its properties will be 
large. At the minimum extreme is a dark disc with pd/ph = 0.25 
and cr rf s lOOkm/s (corresponding to the lowest/highest value 
respectively, found in the simulations 12311 ). Even in this case, 
the annihilation signal from the Earth and Sun are both dom- 
inated by the dark disc rather than the dark halo. However, 
with such a large velocity dispersion, the scaling just described 
means that the dark disc does not lead to the large boosts that 
come from our median dark disc properties. 

The median values for the dark disc properties exclude the 
most probable regions of the CMSSM parameter space. How- 
ever, given the uncertainties in the dark disc properties we can- 
not yet convincingly exclude relevant CMSSM parameter space 
shown in Fig. 4. 

Future surveys of our Galaxy like RAVE [39] and GAIA 
i40ll will detect the local density of dark matter and may disen- 
tangle accreted stars (which will have nearly the same velocity 
dispersion as the dark disc) from those formed in-situ. In this 
case, it will be possible to infer the actual properties of the dark 
disc from these stars, and hence, make more robust predictions 
for the event rate in neutrino telescopes. 

7. Conclusions 

In ACDM, a dark matter disc forms from the accretion of 
satellites in disc galaxies. We have shown how its increased 
phase space density at low velocities enhances the capture rate 
of dark matter particles in the Earth and Sun, resulting in an 
increased neutrino-induced muon flux at the Earth from WIMP 
annihilation. Our main findings are: 

1. The dark disc significantly boosts the capture rate of dark 
matter particles in the Sun and Earth as compared to the 
SHM. This increase owes to the higher phase space den- 
sity at low velocities in the dark disc. For the Sun, the 
expected muon flux from the dark disc with pa/ph = 1 
is increased by one order of magnitude relative to a pure 
SHM-generated flux. If the WIMP is the neutralino in the 
CMSSM, neutrino telescopes will explore a large fraction 
of the CMSSM parameter space. 

2. For the Earth — where WIMP capture and annihilation 
are not in equilibrium — the increase in the muon flux is 
two to three orders of magnitude, although this depends 
sensitively on the distribution function of the dark disc. 
For the SHM alone, the flux from the Sun is far greater 
than that from the Earth. The enhancement from the dark 
disc puts the search for WIMP annihilation in the Earth 
on the same level as the Sun if M x < 100 Ge"vQ. For 
larger WIMP masses, the prospects for detecting muons 



4 Our particular WIMP model already has little parameter space below 100 
GeV, but this extra sensitivity of detection below 100 GeV owes to kinematics 
of Solar system transport and capture, not the WIMP model. 



from annihilation in the Earth requires better models of 
the density of WIMPs bound to the Solar system. 
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